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Abstract 



Under certain mode-matching conditions, small-amplitude waves can be trapped by coupling to soli- 
tons of nonlinear fields. We present a model for this phenomenon, consisting of a linear equation coupled 
to the Korteweg-de Vries (KdV) equation. The model has one parameter, a coupling constant k. For 
one value of the coupling constant, k = 1, the linear equation becomes the linearized KdV equation, for 
which the linear waves can indeed be trapped by solitons, and moreover for which the initial value prob- 
lem for the linear waves has been solved exactly by Sachs [S83] in terms of quadratic forms in the Jost 
eigenfunctions of the associated Schrodinger operator. We consider in detail a different case of weaker 
coupling, K = 1/2. We show that in this case linear waves may again be trapped by solitons, and like the 
stronger coupling case k = 1, the initial value problem for the linear waves can also be solved exactly, 
this time in terms of linear forms in the Jost eigenfunctions. We present a family of exact solutions, and 
we develop the completeness relation for this family of exact solutions, finally giving the solution formula 
for the initial value problem. For k = 1/2, the scattering theory of linear waves trapped by solitons is 
developed. We show that there exists an explicit increasing sequence of bifurcation values of the coupling 
constant, k = 1/2, 1, 5/3, . . ., for which some linear waves may become trapped by solitons. By studying 
a third-order eigenvalue equation, we show that for k < 1/2 all linear waves are scattered by solitons, 
and that for 1/2 < k < 1 as well as for «: > 1 some linear waves are amplified by solitons. 
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1 Introduction 

This paper is concerned with solving the coupled system of equations: 



where k is a real parameter. Of course, the nonlinear equation for A{x^ t) is simply the Korteweg-de Vries 
equation, and it can be solved independently by the inverse-scattering transform [GGKM67]. The coupled 
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system (1) and (2) is a partially linearized version of the system proposed by Hirota and Satsuma [HS81] as 
a model for the dynamics of coupled long waves. 

The coupled system (1) and (2) can be solved exactly when k = 1 and when k= 1/2. The case of k = 1 
is well-known, for then the equation (2) is just the KdV equation itself linearized about the solution A{x, t). 
An elementary exact solution of the linear equation (2) in this case is given by B(x,t) = dxA{x,t). Further 
solutions can be expressed in terms of derivatives of the squared eigenfunctions of the related Schrodinger 
operator with potential A [S83] . 

The case of k = 1/2 is essentially different. In this case, the linear equation (2) is no longer the 
linearization of KdV about any solution. An elementary exact solution of the linear equation in this case is 
given simply by B{x,t) = A{x,t). The main goal of this paper is to construct the general solution of the 
initial value problem for this linear equation when A{x, t) is a multisoliton solution of KdV. 

One way to make clear the difference between the cases k = 1 and k = 1/2 is to consider A{x,t) to be 
the simple soliton solution of KdV (1): 

A{x, t) = I2r]hec\i^{r]{x - ArjH - a)) = -V{x) , (3) 

where x = x — ct — a and the velocity is c = 4?7^. If we look for solutions of the linear equation (2) that are 
traveling waves with speed c, we find the equation 

[-kV{x)B{x)+B"{x)]' = cB'{x). (4) 

Integrating once, using vanishing boundary conditions at x = ±cx;, yields a Schrodinger eigenvalue problem 
for B: 

-B"ix) + kF(x)B(x) = EBix) , (5) 

where E = — c. For k = 1/2, it follows from the fact that B{x,t) = A{x,t) is a solution of (2) that the 
function B{x) = V{x) is an eigenfunction of the Schrodinger operator with eigenvalue E = —c = —'irf. 
Since it has no zeros, it is the ground state eigenfunction. We will see below that there is also one excited 
state for k = 1/2, although it is not relevant here since it corresponds to a different velocity. On the other 
hand, for k = 1, B{x,t) — dxA{x,t) is a solution of (2), which implies that the function B{x) = dxV{x) is 
an eigenfunction of the Schrodinger operator with the same eigenvalue E = —c ~ —Arf. In this case, the 
eigenfunction has a single zero and therefore is the first excited state. It follows that there are at least two 
eigenvalues for k = 1. In fact, there are exactly three states in this case. A final observation is that from 
the construction of the one-soliton solution of KdV (see below) it follows that for k = 1/6, the function 
B{x) — ^(x) is a-n eigenfunction of the Schrodinger operator with eigenvalue E = — c/4 = —if'. It is the 
ground state, and the only eigenfunction. These relationships are summarized in Figure 1. We will have 
more to say about this picture when we discuss the trapping of linear waves by solitons for general values of 
K in § 6. 

The rest of this paper is primarily concerned with developing the general solution of the initial value 
problem for (2) with k = 1/2 when A{x,t) is an A''-soliton solution of KdV (1). In § 2 we show how for 
K = 1/2 a large family of exact solutions of (2) can be obtained from the simultaneous solutions of the 
Lax pair for KdV. When the solution of KdV contains only solitons and no radiation, the construction of 
Lax eigenfunctions is completely algorithmic and algebraic, and consequently the corresponding family of 
solutions of (2) for k = 1/2 can be obtained with great practicality. In § 3 we then establish that in the 
A''-soliton case there are enough of these exact solutions of (2) for n = 1/2 to expand for fixed t any absolutely 
continuous (R) function of x. This fact then leads to a general solution formula for (2) simply by expanding 
the initial data. We present and discuss this formula in § 4. There will turn out to be N linearly independent 
solutions that are asymptotically confined to the union of soliton trajectories and can therefore be considered 
to be bound states. In § 5 we compute the scattering matrix that relates the asymptotic behavior of bound 
states for t — oo to the corresponding behavior for t — + +oo. In § 6 we consider general values of the 
coupling constant k and describe the behavior of some solutions of (2) when A{x, t) is a one-soliton solution 
of KdV (1). These calculations indicate the exceptional nature of the two values n = 1/2 and k= \. In the 
Appendix, we describe several physical applications of the coupled system (1) and (2) to topics in molecular 
dynamics, mechanics, soliton theory, and the fluid dynamics of internal waves. 
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Figure 1: Energy levels of the 12Kr]^sech.^ (rjx) potential for three different values of the coupling constant k. 

2 Exact Solution Formulas for k, = 1/2 

As is well-known [GGKM67] , the KdV equation (1) is the compatibility condition for a pair of linear equations 
involving a complex parameter A for an auxiliary function f{x, t, A). This pair of linear equations is 

Sx/ = -^/-^^/, and dtf = ^d,A-f+(^X^-^A^d,f, (6) 

and is called a Lax pair. A simultaneous solution f{x,t, A) of these linear equations exists if and only if the 
function A{x,t) satisfies KdV (1). Suppose that this is the case. Then, it is a direct matter to verify that 
for fixed but arbitrary A € C, the two functions defined by 



B{x,t) := 



f{x,t,X)e^p ( ±'-{Xx + XH) 



"2 



(7) 



are solutions of the linear equation (2) when k = 1/2. Note here an important point of departure from 
the other solvable case, namely k = 1, where (2) is the linearized KdV equation. In the latter case, 
particular solutions are expressed in terms of the a;-derivative of the square of the Lax eigenfunction f{x, t, A) 
[GGKM74, S83]. By contrast, the formula (7) for solutions of (2) for k = 1/2 is linear in f{x, t, A). This fact 
leads to some important simplifications. 

The formula (7) is only really practical to use if one can explicitly compute the function f{x, t, A). This 
will be the case if the solution A{x, t) of KdV (1) is a pure A''-soliton solution. For each fixed t, A{x, t) is 
then a reflectionless potential of the Schrodinger equation in the Lax pair (6). The multisoliton solutions of 
KdV and the associated solutions of the Lax pair are constructed as follows [KM56]. Let f+{x, t, A) bo given 
by 

/ N-l 

\n-N f / ,\ \ „„„ / I \3+ 



f+ix,t,X) := y-+l. ^^'^ fnix,t)j exp l^--iXx + XH)j , (8) 

where the fn{x, t) are unknown coefficients. Choose N positive numbers rji > r]2 > ■ ■ ■ > rjN , and N arbitrary 
real numbers ai, . . . , ajv, and insist that f+{x, t, A) satisfy the relations 

f+{x,t,2ir,n) = (-l)"+iexp(2ry„a„)/+(x,i,-2ir?„) , (9) 

for all n = 1, . . . , A/'. It is easy to see that these relations imply a square linear algebraic system for the 
coefficients fn{x,t). The determinant of the system is always nonzero, and so the coefficients fn{x,t) are 
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determinc'd uniquely from the sohton eigenvahies {A„ = 2iry„} and the norming constants {a„} in terms of 
exponential functions. From this construction, it can be shown that if one chooses 

A{x,t) :=6id^fN-i{x,t), (10) 

then /_|_(x,t, A) and /_(x,t, A) := — A) arc two simultaneous solutions of the Lax pair (6), and the 

function A{x,t) defined by (10) satisfies KdV (1). The two functions f±{x,t,X) are linearly independent 
for all nonzero A ^ ±2irin. According to the linear relations (9) that determine the coefficients, at the 
exceptional values of A the two functions are proportional. 

The solution A{x, t) of KdV so constructed represents the interaction of A'^ solitons. In particular, as 
t ±oo, the solution can be represented in the form 

N 

A{x,t) ^J2^t{x,t), where At{x,t) := 12r]lsec]i^{r]n{x - a^) - ^vlt) , (11) 

n=l 

where the asymptotic phase constants arc functions of the rjn and «„. 

Below, we will need the asymptotic behavior of the functions f±{x,t,X) as x ^ ±oo for A and t fixed. 
It can be shown that the coefficient functions fn{x, t) remain bounded as x ^ ±oo. Then, letting x tend to 
±00 in the linear relations (9), one finds from dominant balance arguments that 



°° n=0 > 



= 0. (12) 

A=±2ir/„ 



These relations imply that 

CAf-l \ N 

1 + V (±A)"-^/„(ar, t) = A-^ TT (A - 2ir?„) . (13) 
/ 

ra=0 / 71=1 

Therefore, for all A e C, 

lim /+(x,t,A)exp ^(Ax + A3t) =A-^TT(AT2zr,„). (14) 

^ ' n=l 

The large |x| asymptotics for the other solution f-{x,t,X) follow from the definition f-{x,t, A) = f+{x, t,—X). 
From these asymptotics, it is easy to construct the appropriate linear combinations of f±{x,t,X) that corre- 
spond to the Jost functions of the Schrodinger equation, normalized at x = ±oo. 

In the formula (7) we have a choice of sign in the exponent. In fact, it is easy to see that if one considers 
the totality of solutions obtained for all complex A, the choice of sign is redundant. In what follows, we 
adopt a particular choice of the sign and maintain generality by using both Lax eigenfunctions f+{x,t,X) 
and f-{x,t,X). Thus, the particular solutions of the linear equation (2) for k = 1/2 that we will consider 
below will be denoted by h±{x, t, A), given by 



h±{x,t,X) := dxg±{x,t,X) , where g±{x,t, X) := f±{x,t, X) exp ^-{Xx + XH) j . (15) 

Prom the Schrodinger equation (6) for f±{x,t,X) it follows that g±{x,t,X) satisfies the ODE 

-«5xfl± - ^-^9± = >^dx9± ■ (16) 
This ODE plays an important role in suggesting the completeness relation for the solutions h±{x, t, A). 
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3 The Completeness Relation for k = 1/2 



Having in hand a large family of exact solutions of the linear equation (2) for k = 1/2 is certainly useful, but 
wc may then ask whether there arc enough of these solutions to construct the general solution of the initial 
value problem by superposition. A completeness relation is a formula that gives the expansion of arbitrary 
initial data in terms of such a collection of functions. In this section, we will establish the completeness 
relation for the exact solutions h±{x,t,X) obtained in § 2. 

The form of the completeness relation is suggested by a similar argument to that used by Sachs [S83] in 
his investigation of the completeness of squared eigenfunction solutions to the linearized KdV equation. The 
idea is that, ideally we would like to have a differential eigenvalue problem in standard form satisfied by the 
functions h± (x, t, A) : 

L{t)h± {x, t, A) = Xh± {x, t,X), (17) 

where A is the eigenvalue and L{t) is some second-order linear differential operator in x. Then, using the 
two explicit solutions h±{x,t,X) of this problem, we could solve the inhomogeneous problem 

L{t)ip -Xilj = (f) (18) 

by variation of parameters, i.e. by writing as a linear combination of h±{x,t,X) with nonconstant coef- 
ficients, and substituting into (18). For a fixed function (f){x), this determines ip{x,t,X), and we have thus 
constructed the resolvent of the operator L{t), 

ip{x,t,X) = {L{t)-XI)-^(l){x). (19) 

If the spectrum of L{t) is contained in a bounded region of the complex plane (and also under some milder 
conditions), then the Dunford- Taylor integral of the resolvent on a positively oriented contour enclosing the 
spectrum yields the identity operator: 

--^ ii}{x,t,X)dX = -^. (f{L{t)-XI)-'^(t){x)dX = <j){x). (20) 

znl J ZTTl J 

But we do not have a second-order eigenvalue problem for h±{x,t, A). Instead we have the second-order 

equation (16) for g±{x,t, A). However, we make the guess that a similar procedure will apply here. Namely, 
for appropriate side conditions (see below) we solve the inhomogeneous equation 

-ia^^-i^V-Aa,V = </', (21) 

for ?/;(a;, t, A) using variation of parameters with the two functions g±{x, t, A) solving the homogeneous equa- 
tion (16), and then we differentiate the resulting formula with respect to x. Formally speaking only, we have 
thus constructed the resolvent of the "operator" 

L{t) = -id, - i^d-^ . (22) 

The obstruction to rigor here is that is not well-defined. Nonetheless, we are guided to hypothesize that 

-^.jd^i^ix, t, A) dX = (l>{x) , (23) 

for an appropriate contour of integration. This formula turns out to be correct, although a direct proof must 
be supplied. The proof we use follows Miller and Akhmediev [MA98]. 

3.1 Solving the inhomogeneous problem. 

We express the solution of the inhomogeneous problem in the form 

i){x, t, A) = C+ {x, t, X)g+ {x, t, A) + C_ (a;, t, X)g- {x, t, A) , (24) 
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subject to the usual "reduction of order" condition 

d^C+ {x, t, A) -g+ix, t, A) + d^C_ {x, t, X) ■ g_{x,t,X) = . (25) 

Substituting (24) into the equation for tp, and using (25), one finds 

a r< f + \^ ■Hx)g-{x,t, X) a n i ^ w ■(t>{x)9+{x,t,X) 

d:,C+{x,t, A) = -«^Tj77 , and a^C_(a;,t, A) = i^tjtj ^ , (26) 

where W{g+, g-) := g+dxg- — g-dxg+ is the Wronskian. 

From the differential equation (16) satisfied by g±{x,t, A), it follows that 

dxW{g+,g^)=iXW{g+,g_). (27) 

Using the large asymptotics of f±{x,t,X) obtained in § 2 one then solves (27) uniquely and finds that 

N 

W{g+, g.) = iAi-2^ exp(z(Ax + XH)) [] (A^ + Ar^) . (28) 

n=l 

In solving the inhomogeneous equation (21) for t/', we really should impose appropriate side conditions. Here, 
the side conditions we use are not related to boundary conditions in x as much as analyticity conditions in 
A. It is easy to check that for each a;o,c/, the function ipu{x, t, X) defined by 



n=l 



+ r ^.(^.^.A)expH(A. + A3t))^(.)^^.^_^^^^^^^^ 



(29) 



is a solution analytic in A for 9(A) > and |A| sufficiently large. Similarly, il)L{x,t,X) defined for each xo,l 

(30) 

,+ (., t, A) exp(-i(A._+A3tM£) . ^ 



f 

J X 



N 



X'-"^l[{X'+4rjl) 



is a solution anal}d;ic for Sj(A) < and |A| sufficiently large. The qualification of |A| being sufficiently large 
is necessary because the expressions have poles at the soliton eigenvalues in the respective half-planes where 
the two functions g±{x,t, A) become proportional. However, these are the only finite singularities and both 
solutions ipu{x, t, A) and iI)l{x, t, A) are meromorphic in the whole of their respective open half-planes. 

The arbitrariness of the parameters xo,u and Xo,l would seem to be a problem; but it will turn out 
that these terms contribute nothing to the Dunford- Taylor integral that we will prove gives the required 
completeness relation. 

3.2 Integrating the resolvent. 

Here we show that the guess we made is indeed correct. 
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Theorem 1 Let (f){x) be an absolutely continuous function in L^(M). Let xo,u o,nd xq^l be constants and let 
f G M 6e fixed. Then, 



(b(x) = lim 



dxtlJu{x,t,X)dX+ / dxtpL{x,t,X)dX 



(31) 



where Cu is the positively oriented half-circle from R to —R in the upper half-plane, and Cl is the positively 
oriented half-circle from —R to R in the lower half-plane. 

Proof: First, we show that the terms depending on the arbitrary parameters xo,u and xq^l converge to 
zero as R^ oo. This will justify calling the function ipu or ipL a "resolvent" even though the inverse is not 
unique. Consider the integral 



Ju ■= / da: 
ICu 



I 



g-{z,t,X)exp{-i{Xz-\-XH))(f){z) 



N 



dz-g+{x,t,X) 



X^-'^l[{X^ + A,l) 



n=l 



dX 



(32) 



Ici 



X'''h+{x,t,X) 

N 



Q- {z, t, A) exp{—i{Xz + X^t))(l){z) dz dX , 



n, o ■J3:o,u 
(A' +477^) 

n=l 

where we have used the relation (25). Recall that 

h+{x,t,X) = ^X~'^dxfN-n{x,t) , 



N 



n=l 



(33) 



TV 



Tl=l 



<?_(z,i,A)exp(-i(Az + A3i)) = i + ^(-A)-'7Ar_„(z,i). 
Therefore, for all A with |A| = i? > 1, 

1 ^ 

\h+{x,t,X)\ <-^'Yl \dxfN-n{x,t)\ , 

and 



N 



su-p\g_{z,t,X)exp{-i{Xz + X^t))\ < 1 + sup ^ |/jv-„(2;, i)| 



(34) 



(35) 



n=l 



This latter relation assumes the uniform boundedness of the functions fk{z,t) in z. Finally, it is clear that 
for |A| = i? > sup„2?7„, 

Aril 



It follows that for all A with |A| = i? sufficiently large 



N 



>i?n(i 



n=l 



i?2 



\Ju\<'^Uh, 



where 



N 



71=1 



K{x,t)=lTX{ (l-^) ■[Y.\d,fN-n{x,t)\] ■ (l + SUp^|/^_„(z,i)| 



-1 / N 



N 



(36) 



(37) 



(38) 
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The bound (37) clearly vanishes as i? — > oo. A nearly identical argument shows that the integral 



Jl 



Q- {z, t, A) exp{-i{\z + X^t))<p{z) 



N 



dz-g+{x,t, A) 



n=l 



dX 



(39) 



satisfies the same bound (37) as Ju- 

Now we consider integrating the second terms of dxtpu{x, t, A) and dxi)L{x, t, A) respectively. For brevity, 
define 

g+{z, t, A) exp(— z(A^ + X^t))h-{x, t, A) 



Yix,z,t,X) :-- 



N 



(40) 



X'-'^^HiX'+Arjl) 



n=l 



Note that this can be written as 



Y{x, z, t, A) = exp(iA(a; — z))- 



1 + E 

^ n—l 



fN-n{z,t) 

A" 



An 1+ 



n=l 



A2 



(41) 



^x[l+Y: 



N 

E 



dx!N-n{x,t) 

(-A)" 



and therefore, 



Y{x, z, t,X) = i exp{iX{x — z)) (1 + A(a;, z, t, A))) , 



(42) 



where A(x, A) ~ 0(A^^) uniformly in x and z for fixed t. It also follows from additional cancellation 
that for z = X, A{x,x,t, X) = 0{X~^) uniformly in x. Finally, derivatives of A are controlled as well: 
dzA{x, z, t, A) — 0{X^^) uniformly. The integral we need to compute for the contribution of dxi^u{x, t, A) is 



/ / Y{x,z,t,X)(j){z)dzdX = i / / exp{iX{x — z))(j){z) dz dX 

JCh''—^ JCrJ—oo 

+ i / A{x,z,t,X)exp{iX{x — z))(l){z)dzdX. 

Cf!. J —OQ 

Note that since the integrand is analytic in the upper half-plane, the first term can be written as: 



(43) 



exp(iA(a; — z))(l){z) dzdX = —i 



Cr J —oo 



R px 



exp{iX{x — z))(l){z) dz dX . 



(44) 



B.J-00 



In order to control the error term, it is necessary to integrate by parts once: 



i / A{x,z,t,X)exp{iX{x — z))(l){z)dzdX = —(j){x) I 
exp(iA(a; - z)) 



A(x, X, t, A) 



dX 



Cr 



(45) 



-a^(A(x, z, t, X)(t){z)) dz dX . 



Cr J -oo 



The boundary term at z = —oo vanishes because A{x, z, t, A) is bounded there, (j) is continuous and integrable, 
and exp(iA(a; — z)) is exponentially small for SJ(A) > 0. Since the exponential is bounded in magnitude by 
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unity for 5(A) > 0, and the contour is of length ttR, the above estimates of A imply that there exist Ko{x, t), 
Ki{x,t) and K2{x,t) all positive, such that 



A(a;, z, t, A) exp(zA(a; — z))(t){z) dz dX 



Cr J -oo 

Ko{x,t) 



Ki{x,t) ^ 
R 



(46) 



This proves that 



lim 

iJ->CXD 



Cr J -co 



Y{x, z,t,\)(j){z) dz d\ = —i Um 



R — >oo 



exp(iA(a; — z))(f}{z) dz dX . 



RJ -co 



Similar arguments applied to the contribution of dxipL{x,t, X) show that 

rR POO 



r rco rri pco 

— lim / / Y{x,z,t,X)(l){z)dzdX = —i hm / / exp{iX{x — z))(l){z) dz dX , 

Jcl Jx J-R Jx 



and therefore 



lim 

2m R^oo 



Cr 



dor:'il^u{x,t,X)dX+ dx'ipL{x,t,X)dX 



Cl 



(47) 



(48) 



(49) 



1 



277 R- 



lim 



exp(iA(x — z))(j){z) dz dX = 4){x) , 

-R J-oo 

with the last equality following from Fourier inversion. This establishes (31) and the theorem. □ 

As it stands, the completeness relation given in Theorem 1 is not really an expansion of (j){x) in terms of 
the functions h_{x,t, A), because the expansion coefficients themselves depend on x. This is easily remedied 
by casting the right-hand side of the completeness relation into a more useful form. We do this now. 

Theorem 2 Let 4>{x) be an absolutely continuous function in L^(R). Let <; G M &e fixed, and choose any 
t« G M U {— oo, +oo}. Define the mode function 

H{x, t, A) := X^h- {x, t, A) , (50) 

which is an entire function of X, and the amplitudes 

X^g+{z, t, A) exp(-i(A^ + X^t)) 



b+{t,X) := / 



N 



b-{t,X) :-- 



A n (A' +477^) 

n=l 

^ X^ g+{z,t,\)cyi^{-i{Xz + XH)) 



CO 



N 



{z) dz , 



(j){z) dz , 



(51) 



xJlix^ + 'inl) 



71=1 



and set b{t,X) := 6+(t, A) + b {t,X). The amplitudes have simple poles at X = and X = ±2ir]n for 
n = 1, . . . ,N . Finally, set 



&o(t) := i Res (6+(t,A) -6"(i, A)), and b^{t):=T Res b'^{t,X). 

A=0 X=±2iri„ 



(52) 



Then we have the expansion 



1 



R^oo 2wi 



0(x) = lim — P.V. / b{t,X)H{x,t,X)dX 



N 



(53) 



+ bo{t)H{x,t,Q) + J2[Kit)H{x,t,-2ijjn)+b+{t)H{x,t,2irin)] ■ 



n=l 
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Remetrk: Since w is now fixed and not a function of x, this expansion (53) is a true completeness relation, 
expressing an arbitrary given function (p{x) as a sum of known functions H{x, t, A). From the exact formulas 
(50), (15), and (8), it is clear that the part of the expansion (53) represented by the singular integral is Fourier- 
like, with the corresponding components of the solution, H{x,t,X) for A S K being bounded oscillatory 
functions tending to complex exponentials for large x. On the other hand, the discrete contributions to the 
solution represent bound states. The 2A'' + 1 bound state terms in (53) are not linearly independent. From 
the fact that at the eigenvalues ±2ir]n the functions g- (x, t, A) are all linear combinations of the same A'' 
functions fo{x, t), . . . , fN-i{x, t), it is clear that only N of the bound states are linearly independent. These 
facts are easiest to see when one takes w to oo or — oo. Then, half of the contributions from the eigenvalues 
disappear, and it only remains to express the bomid state at zero, H{x, t, 0), in terms of H{x, t, ±2irin)- This 
can be done directly. From the exact formulas (50), (15) and (8), we see that 

H{x,t,0) = {-lfd,foix,t), (54) 

and making use of the relations (9) satisfied by /+ at the eigenvalues, 

N-l 

H{x, t, 2ir]n) = (-1)"+^ exp(-2?7„a„) ^ {2ir]nTdJp{x, t) . (55) 

Expressing H{x, t, 0) in terms of H{x, t, 2ir]n) is therefore a polynomial interpolation problem. Introduce the 
polynomial 

N-l 

P{X) = d.fp{x,t)XP . (56) 

p=0 

Given isolated values of this polynomial: 

P(2i77„) = (-l)"+iexp(2r/„a„)F(x,t,2ir;„), for n = l,...,N, (57) 

we are to find P(0) and thus H{x, t, 0) = {—l)^dxfo{x, t) = (— l)^P(O). Expressing P(A) explicitly in terms 
of Lagrange polynomials gives 

PW = E(-l)"^' eM'^Vnan)H{x, t, 2i7?„) [] ^ , (58) 

n=l k^n 



and therefore 



N 



Hix,t,0) = J2 



(-l)"exp(2?7„a„) [] 



Vk 



Vn - Vk 



H{x, t, 2ir]n) . (59) 



There is of course a similar expression for H{x, t, 0) in terms of H{x, t, —2ir]n)- □ 

Remark: An important distinction between the completeness relation stated in Theorem 2 and that found 

by Sachs [S83] for derivatives of squared Schrodingcr cigcnfunctions is in the nature of the singularity at 
A = 0. Sachs shows that in the expansion of (p in terms of derivatives of squared eigenfunctions, there is an 
apparent singularity at A = that it is in fact removable. On the other hand, the integral in Theorem 2 is 

essentially singular and the residue contribution of bo{t) is nonzero. □ 

Proof of Theorem 2: We first establish that in the formulas (29) for dxi^uix, t, A) and (30) for dxi^hix, t, A) 
we may replace x in the limits of integration by any other value without changing the result of the theorem. 
That is, we will now show that the integral 

Y{x,z,t,X)(p{z)dzd\- I I Y{x,z,t,X)(j){z)dzd\, (60) 

Cjj J — OQ JCl 
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which wc have ah'cady seen converges as R tends to infinity to —2'Ki(p{x) in the ease that w = x, is in fact 
independent of w. Holding R fixed and differentiating with respect to w, we must show that for sufficiently 
large R, 

(j){w) (L Y{x,w,t,\)d\ = Q, (61) 

J\\\=R 

identically in x, w, and t. Being as the integrand is meromorphic in the finite A plane, we can evaluate the 
integral by residues. There are simple poles at A = and A = ±2irin (ot n = 1, . . . ,N. Using the linear 
relations (9) satisfied by f± at the eigenvalues A = 2jry„, we find 

Res Y{x,w,t,X) = Y.^^^ d,fp{x,t) 

N-1 



p,9=0 

(62) 

Res Y{x,w,t,X) = ^(-l)^-fii^^i- d.U{x,t) 



p,q=0 " 



where 



N 

Dk ■= n (2«% - 2jr?n) n + ^^J?") • (63) 

Similarly, for the residue at zero. 



Res y(x,«,,t,A) = (-l)^Mf!l|MM. (64) 

n 



Adding all the residues and collecting coefficients of the terms dxfp{x, t) and fg{w, t)dxfp{x, t), we find that 
the sum of the residues will be zero if: 

^ (2i?7fc)P 

/p.= ^-^=0 (65) 
k=i " 

for all odd p= 1, 3, 5, ... , 2A'' — 3, and if: 

1 ^ 1 

-Tr + ^=0- (66) 

n=l 

These expressions are themselves sums of residues of meromorphic differentials. Thus, by inspection, one 
finds that for p = 1,3,5,..., 2N — 3, 

1 f XPdX 
If 



n=l 
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where C is any simple c;ounter-clockwise oriented contour that encircles the points A = 2irjn for n = 1, . . . ,N 
(but without enclosing the conjugate eigenvalues or A = 0). With p bounded by 2N — 3, the path of 
integration can be blown out to infinity in the upper half-plane, and then brought down to the real axis so 
that 

1 /-^ \Pd\ 
= ^''^ 

with the last equality following from the oddness of the integrand for odd p. Finally, consider the integral 
I -I defined by 

n=l 

Evaluating the residues inside C, we find 



^- = E • (70) 



fc=i 



On the other hand, we can again blow the contour C out to infinity in the upper half-plane and bring it 
down to the real axis. This time, there is a singularity at A = 0, so the Plemelj formula must be used. We 



Once again, by oddness, the principal value integral vanishes identically, and then combining this result with 
the previous expression we obtain the required vanishing. 
This shows that for any w, 



p PW p PCO 

/ / Y{x,z,t,X)(t){z)dzdX~ / / Y{x,z,t,X)(l){z)dzdX 

'J Cu — CO J Cl '^w 



(72) 



Establishing (53) and therefore the Theorem now amounts to using the residue theorem once again to 
deform the integration paths Cu and Cl in (72) to the real axis. One finds discrete contributions at 
the poles A = ±2i77„ and then applying the Plemelj formula to contract the; contour to the real axis in 
the neighborhood of A = gives a discrete contribution proportional to H{x,t,0) and the principal value 
regularization of the singular integral over the continuous spectrum. □ 



4 Solution of the Initial Value Problem for k, = 1/2 

It is easy to see that when A{x, t) is an 7V-soliton solution of KdV (1), one can use the completeness relation 
to solve the initial value problem: 



dtB + 



AB + dlB 



0, 



B{x,Q) = (t){x). 



(73) 



Setting t = Q, and picking a convenient value of w, say, w = +oo, one computes the amplitudes (51) and 
discrete coefficients (52). Then, because the function H{x,t, A) satisfies (2) for k = 1/2 and for each complex 
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A, the expression 



B(x,t) := lim — P. V. 



b{0,X)H{x,t,X) dX 



R 

N 



(74) 



provides the solution of the initial value problem (73), generally in the sense of distributions. That is, 
B{x, 0) = 4){x) by Theorem 2, and for each test function ^p{x, t) that is differentiable in t and three times 
differentiable in x and has compact support in (a;,f) G M x M+, one shows by exchanging the order of 
integration that 



oo POO 



dtip{x, t) + - A(x, t)d^ip{x, t) + dlip{x, t) 



B{x, t)dxdt = 0. 



(75) 



The solution will be classical in as much as it is possible to differentiate with respect to x and t under the 
integral sign in the solution formula (74). This requires additional smoothness and decay assumptions on 
the initial data (l){x) that we do not consider here. 



5 Scattering of Bound States for k, = 1/2 

Of particular interest in applications is the A^-dimcnsional (recall that A{x, t) is an A'^-soliton solution of 
KdV) subspace of solutions of (2) for k = 1/2 consisting of bound states. This subspace represents linear 
waves that are trapped by the solitons of the potential A{x, t). For large \t\, these bound state solutions are 
all confined to the trajectories of the solitons. Therefore, it follows that each bound state B{x, t) has two 
asymptotic representations: 

N 

B{x,t)^Y.PnAt{x,t), i^ioo, (76) 

n=l 

for some constants depending on B{x,t), where A^{x,t) are defined by (11). Since there are exactly 
A'^ linearly independent bound states, it follows that the constants are completely detc^rniincd from the 
constants /?~. In particular, there exists an invertible N x N matrix T with entries depending only on the 
data specifying the A/'-soliton solution A{x,t), such that 

JV 

fe=i 

The matrix T is called the bound state scattering m,a,trix. In this section, we compute the scattering matrix 
explicitly, and show that its elements only depend on the soliton eigenvalues 771, ... , 77 at. 

If A{x, t) is an A^-soliton solution of KdV (1), then a family of solutions to (2) for k = 1/2, parametrized 
by complex A, is given by 

A(x t) 

h+{x, t, A) = -gl^ + ^ X^'-^dJnix, t) . (78) 

n— 

We want to analyze these solutions in the limit of large \t\, in a frame of reference traveling with constant 
velocity c. 

The first step is to see how the coefficients /„(.t, t) behave for large \t\. Let x = .t — ct be fixed as t = t 
goes to either +00 or —00. Begin by taking 77^ < 4c < ?7^-i to see how the coefficients behave in between 
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the solitons. In the limit r — > +00, the equations (9) imply that 



1 + ^(-2z,7„)'=-^/fe 



0, n = 1, . . . ,m — 1 . 



k=0 
N-l 



(79) 



1 + 5^(2ir?„)'=-^/fc ^ 0, n = m,...,N. 



This is an invertible Vandermonde system for the coefScients fk, so that as r +00, the fk all become 
constants, independent of % and r. Thus, dxfk{x,t) vanishes between the solitons for all k. The analogous 
result holds as t — 00. This shows that the solutions of (2) for k — 1/2 described by the formula (78) 
are asymptotically confined to the individual frames of reference of the moving solitons in the potential field 
A{x,t). 

Now set c = 4?7^ to go into the moving frame of reference of one of the solitons. Taking the limit r — > +00 
yields 



Af-l 



1 + ^(-2ir?„)'=-^/fc 



0, n = l,...,m— 1, 



fc=0 
JV-l 



(80) 



1 + ^ (2ir?j'=-^/fc ^ 0, n = m + l,...,N. 



k=0 



This is a system of — 1 equations in A'^ unknowns, so it can be used to asymptotically eliminate d-^fo 
through in favor of O-^/n-i, which we know is proportional to the A''-soliton solution of KdV, A{x, t). 

Thus, as T — > +00, with c = 4?7^, 



dyfk — Qmk9-yfN-l — 'Tr^QmkA , 

6i 



(81) 



for k = 0, . . . , N — 2 where the numbers Qmk are the unique solution of the inhomogeneous system of linear 
algebraic equations 



N-2 



(-2ir?„)-i + ^ (-2ir?„)'=-^Q„fc = 0, n=l,...,m-l, 



fe=0 
N-2 



(82) 



(2ir?„)-i + ^(2i77„)'=-^Q„fe = 0, n = m + l,...,N. 



fe=0 



One can similarly show that as r ^ — 00, with c = 4?7^, 



(83) 



for /c = 0, . . . , iV — 2, where the star denotes complex conjugation. 

Now consider particular solutions Bj{x,t) of (2) for /t = 1/2 obtained as linear combinations of A'' others 
expressed by the formula (78) evaluated on the N soliton eigenvalues. The formula for Bj{x,t) is 



N 



N 



Bj{x,t) = ^Fjkh+{x,t,2ir]k) =^Fjk 



k=l 



k=l 



A{x,t) 



N-2 



+ ^(2%)"-^5,/„(:r,i) 



n=0 



(84) 



with c = 4//^ 



where F = {Fjk} is a matrix of arbitrary constants. From the asymptotics of fnix, t), we have as r — > —00 
^ A > ' Fii^GZ^ , where GZ^ := + ^ > ' (2wfc)"-™Ql„ • (85) 



^Y^F^kG-,^^ where := + 1 ^^(2z%)"-^Q: 



k=l 
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So, with the choice that the matrix {Fjk} is the inverse of the matrix G = {G^^^^}, the particular sohition 
Bj{x, t) of (2) for K = 1/2 will be completely confined as t — > — oo to the frame of reference moving with speed 
c = irjj, where it will be locally indistinguishable from the solution A{x,t) of KdV. Let us now determine 
how Bj{x, t) will behave in the various soliton frames as t +00. Passing to the limit of r — > +00 in the 
frame with velocity c = 4r]^ gives 

JV N-2 

Bj^AY,F,,G+^, where G+„ :=--— + - J] (2z,yfc)"-^g„n • (86) 

These asymptotics give us a formula for the bound state scattering matrix: 

T:= [(G")-iG+]^ . (87) 

It is clear that the elements of T depend only on the N soliton eigenvalues r]i, . . . ,r]N- There is no dependence 
on the soliton phase variables ai, . . . ,aN. Therefore, the asymptotic scattering properties of linear waves in 

(2) with K = 1/2 are insensitive to phase shifts among the solitons in the potential A{x,t). As a concrete 
example of the scattering matrix, we compute it explicitly for = 2 for arbitrary ryi > 772 > 0: 



? 2 



(88) 



The fact that T22 is negative means that it is possible for the interactions of the solitons in A(x, t) to convert 
trapped linear waves of elevation into waves of depression, and vice-versa. 



6 General Values of k 

We expect that for most values of k, the linear waves satisfying (2) will not be permanently trapped by solitons 
present in the potential A{x, t). This is suggested by considering the simplest case, namely taking A{x, t) to 
be the one-soliton solution of KdV (1). The soliton travels with velocity c = 477^, so that A = —V{x) with 
X = X — ct — a. Corresponding traveling wave solutions B{x) of the linear problem that propagate with the 
same velocity and decay as x — > ±00 satisfy 

-B"{x) + i^V{x)B{x) = -cB{x) . (89) 

Since c is fixed, we can view this as an eigenvalue equation with k as the eigenvalue. We therefore expect that 
only isolated values of k will admit nontrivial decaying solutions B{x)- We have already seen that k = 1/2 
and K = 1 are indeed eigenvalues. For k = 1/2 the eigenfunction B{x) is an even function of Xi while for 
K=l the eigenfunction B{x) is odd in x- Since eigenfunctions of (89) must be nondegenerate and therefore 
have either odd or even parity in x, there cannot exist a nontrivial bound state eigenfunction of (89) for all 
«; e [1/2, 1] since the eigenfunction would have to change parity from one endpoint to the other. Therefore, 
at least one value of k e [1/2, 1] is not an eigenvalue. For such k, there is no bound state traveling wave 
solution of (2) that is trapped in the soliton trajectory. 

We can be more precise about this phenomenon. The left-hand side of (89) can also be viewed as a 
Schrodinger operator L{k) depending on a coupling constant k, and the condition for wave trapping by 
solitons is simply that — c £ Sp(L(K)), where Sp denotes the point spectrum. The number of discrete 
eigenvalues is a nondecreasing function of k > 0, corresponding to the deepening of the potential well. There 
is an infinite unbounded sequence of cutojf\ahies k^"* of k at which the number of eigenvalues changes by 
one, and the new eigenvalue is born from the continuum. Each eigenvalue, once born, is distinct and is a 
decreasing function of k. From these arguments, it follows that there exists an infinite unbounded sequence 
of bifurcation values of k at which one eigenvalue crosses the level E = —c, and a bound state traveling 
wave solution of (2) exists. 

It is easy to find the bifurcation points because the hyperbolic secant squared potential is so well under- 
stood. The potential kV{x) is exactly reflectionless for 12k = 71(71+1) for n= 1, 2, 3, The corresponding 

energy levels are En,k = —k'^rf, for A; = 1, . . . ,n. Therefore, for n > 1 in this sequence, there is always one 
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eigcnvahic that is exactly equal to — c = —At]^. The corresponcling eigenstate is always the n — 1st state 
and therefore has n — 2 zeros. It follows that the bifurcation points are k = kJJ'' = (n + l)(n + 2)/12 for 
n= 1,2,3,.... 

The fact that some linear waves may be permanently trapped by isolated solitons at a bifurcation point 
K = does not necessarily imply that there will be no losses to radiation when solitons in the field A{x, t) 
interact with one another. Such a lossless interaction might suggest the "integrability" of the linear equation 
(2). We have indeed seen that this is the case for the first two bifurcation points, k = n\^^ = 1/2 and 
K = ^2'^ = 1, but it is by no means clear that the trend continues for higher-order bifurcation points. For 
the rest of this section, we therefore restrict attention to the case N = 1, that is, we take the nonlinear field 
A{x,t) to be a one-soliton solution of KdV (1). 

Using (3) and the change of variables x' = r]{x — 477^f — a) and t' = rft, the equation (2) becomes, after 
dropping primes, 

dtB + d^[-AB + 12«;sech^(a;)S + dlB] = . (90) 

This equation is of course solved by separation of variables. We seek separated solutions B{x, t) = ba{x) exp(crt) 
and obtain the third-order eigenvalue problem 

[Ab^ix) - 12Ksech^{x)b^{x) - b'^{x)]' = (jb^{x) , (91) 

where the prime denotes differentiation with respect to x. In this context, what we have been calling "trapped 
linear waves" correspond to bound-state eigenfunctions of (91) with a ~ 0. Such solutions have finite mass 
and energy and are stationary in the moving frame of reference of the soliton A{x,t). As we know, such 
eigenfunctions with (7 = exist only at the bifurcation values of k = kJJ'^. However, it is clear that for general 
values of k there are other possibilities. There may be eigenvalues a that are purely imaginary, giving rise 
to oscillating modes that travel in the soliton frame. More generally, if an eigenvalue has a nonzero real part 
for some k, then there will be a mode that is either amplified or exponentially damped as it propagates with 
the soliton. 

The eigenvalue problem (91) has two simple symmetries. Whenever ba-{x) is an eigenfunction with 
eigenvalue a, then bcr{—x) is an eigenfunction with eigenvalue —a and ba{x)* is an eigenfunction with 
eigenvalue a* . Therefore, the eigenvalues either come in purely real pairs (|cr|, — |ct|), purely imaginary pairs 
(i|cr|, — i|(T|), or in complex quartets (a, — a, cr*, — cr*). These symmetries indicate the distinguished role of 
cr = as a point that if it appears in the spectrum for some k can signal a bifurcation in the number of 
eigenvalues. This explains our terminology and notation for the values k = . 

Most points on the imaginary cr axis correspond to continuous spectrum. This can be seen by the following 
argument. Let k be fixed. Suppose a = iuj with a; e M. For large |a;|, the solutions of (91) have the form 
of linear combinations of eyi\>{ikijjjx) where k = k^ij are the three roots oi k^ + ik — w = 0. Exactly one of 
these roots, say k = ki^^, is real, while the other two form a complex-conjugate pair. If we seek a generalized 
eigenfunction normalized to exp{iku,flx) as a; — » —00 through a "shooting" method, we have three complex 
constants to exploit: the coefficient of the decaying mode for large negative x, and the coefficients of the 
decaying mode for large positive x and the finite amplitude contribution for large positive x. Matching the 
values of bia,{x), b^^{x), and b'l^{x) at a; = gives three complex equations in three complex unknowns. 
If this system of equations is solvable at all, one expects it to be solvable for almost all real u>, yielding a 
generalized eigenfunction. For exceptional values of cj where there is not a generalized eigenfunction, there 
will be a genuine bound-state eigenfunction, since the spectrum is a closed set. 

We have used a numerical Fourier-based collocation (pseudospectral) method to find the discrete eigen- 
values of (91) over a range of values of the coupling constant k. Essentially this involves approximating 
the continuous function bcr{x) by a periodic discrete series. Then the derivative dx can be approximated to 
exponential accuracy by a derivative matrix D, for which an explicit formula is given in [CHQZ88]. This 
is then used to construct a discrete approximation to the operator on the left-hand side of (91). Standard 
techniques can then be used to obtain the eigenvalues and eigenvectors of this matrix. The corresponding 
eigenfunctions always decay exponentially, but sometimes they decay very slowly for large x of one or the 
other sign — luckily not both. To obtain accurate results it was necessary in these cases to change the 
dependent variable by multiplying by an appropriate exponential function of x to enhance the decay on 
the slowly-decaying side without changing decay into growth on the other side. Our results over the range 
< «; < 5 are shown in Figure 2. The bifurcation values k^'^ appear to be of two different types. If n is 
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Figure 2: Real (above) and imaginary (below) parts of the discrete eigenvalues a for the eigenvalue problem 
(91) as a function of the parameter k. Different eigenvalue branches are displayed with different styles of 
lines (solid, dashed, etc.). 



odd. then when k increases through the value k = k^^^ , a new pair of real eigenvalues is born from the origin 
a = 0. As K is further increased, the pair of eigenvalues moves at first away from the imaginary axis and 
then changes direction and contracts toward the origin. When k increases through the even bifurcation value 
K = 't^+i, the pair enters the origin and re-emerges as a complex eigenvalue quartet. Further increasing the 
value of K causes the quartet of eigenvalues to move through a maximum in the magnitude of the real part 
and then toward the imaginary axis with the magnitude of the real part decreasing to zero while the magni- 
tude of the imaginary part increases without bound. It does not appear that the quartet of eigenvalues ever 
disappears into the continuous spectrum, although it comes arbitrarily close as k increases. This scenario is 
repeated again and again as k increases through each odd bifurcation value. Representative eigenfunctions 
are plotted in Figure 3. Here, one can see that when the eigenvalue a has a nonzero imaginary part, the 
decay of the eigenfunction can be quite slow on the "downstream" side of the soliton. As remarked above, 
this effect is compensated for in our numerics by working with a modified eigenfunction. 

There are no discrete eigenvalues at all for k < 1/2 (and in particular for /t < 0), and for all k satisfying 
1/2 < K < 1 and k > 1, there is always at least one eigenvalue with a nonzero real part, which corresponds 
to an exponentially growing eigenfunction and therefore instability. The values k = 1/2 and k = 1 are 
distinguished as the only values for which there exist discrete eigenvalues and at the same time all eigenvalues 
have zero real parts, so the system is neutrally stable. For all other values of k, either there are no discrete 
eigenvalues at all in which case all initial conditions for (2) disperse away algebraically in time, or there are 
discrete eigenvalues with positive real parts in which case the linear waves are amplified by the soliton in 
the field A{x,t). 
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Figure 3: Left: the real-valued eigenfunction corresponding to the eigenvalue a with positive real part for 

K = 0.85. Right: the com,plex-valued eigenfunction corresponding to the eigenvalue a in the first quadrant 
for K = 1.2. The solid curve is the real part and the dashed curve is the imaginary part. 

7 Conclusion 

The coupled system consisting of the KdV equation (1) and the linear equation (2) is integrable for two 
distinct values of the coupling parameter k. The integrable case of k = 1 has been studied by other authors 
[GGKM74, S83]. In this paper, wc have given new results for the other integrable case, namely n = 1/2. 
In particular, we have shown how to construct the general solution of (2) for k = 1/2 when the nonlinear 
field A{x,t) is an A''-soliton solution of the KdV equation. This general solution is represented in terms of a 
number of bound states (equal to the number N of solitons in the field A{x, t)) and a contimmm superposition 
of radiative states given by a singular integral. With the help of numerical computations, we have placed the 
integrable cases in context by examining the behavior of the linear equation (2) for general values of k, when 
A{x,t) is a one-soliton solution of KdV. These calculations show that the linear equation (2) behaves as an 
unstable dynamical system for most positive k. The integrable value of k = 1, for which the equation (2) is 
the linearized KdV equation, is an isolated stable point, since a small change of either sign in the value of k 
will lead to the presence of exponentially growing modes. The other integrable value of k = 1/2 represents 
the boundary between a stable system without any bound states for k < 1/2 and an exponentially unstable 
system for k > 1/2. 

In physical applications of the coupled system (1) and (2) as discussed in the Appendix, the presence 
of instabilities indicates that more terms need to be included in the model. However, in the stable cases 
the model is indeed expected to be physically meaningful. And in this regard, the two integrable cases can 
provide useful starting points for perturbation theory. 

As a final remark, let us indicate the kind of calculations that are possible for the coupled system (1) and 
(2) for K = 1/2 with the aid of the completeness relation. For a family of relevant initial data for the linear 
equation, one can explicitly compute the projection onto the bound states, and consequently determine the 
long time behavior of the corresponding solution of (2). Also, the long time behavior of the dispersive part of 
the solution can be computed from the explicit representation of this component of the solution as a singular 
integral. We leave such applications of the completeness relation for further investigations. 

Appendix: Some Applications 

It is useful to keep in mind some applications in which the coupled system (1) and (2) might arise. In fact, such 
equations appear in the modeling of coupling of acoustic phonons in long polymer molecules. Many organic 
polymers {e.g. DNA and a- helix proteins like acetanelide) may be considered from the mechanical point 
of view as long chains of nearly identical masses. This "backbone" of the molec\.ile supports a longitudinal 
vibrational mode in which the masses are all moving in tandem with zero frequency {i.e. simple translation) 
in the long-wave limit; the associated quanta are called acoustic phonons. In the presence of intrinsic weak 
nonlinearity, the Korteweg-de Vries equation describes these vibrations in the long- wave limit. Usually, the 
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masses making up the chain contain internal degrees of freedom (e.g. the "breathing" modes of base-pairs 
in DNA, and the so-called amide I exciton modes of the C=0 bond in each peptide group of an a-helix 
protein). The coupling of these internal degrees of freedom to the motion of the backbone leads to a variety 
of interesting dynamical models {e.g. the discrete sine-Gordon equation for DNA, and the discrete nonlinear 
Schrodinger equation for a-helix proteins). 

We may consider a situation in which the internal degrees of freedom are themselves acoustic phonons 
associated with transverse vibrational modes. This can be visualized with the help of a concrete mechanical 
model, whose equilibrium configuration is shown in Figure 4. The backbone is made of heavy masses M 




L ^ M 



Figure 4: The equilibrium configuration of the mechanical model. 

connected by stiff springs. Mounted on each heavy mass is a transversely-oriented frictionless track in which 
rides a small mass m. The mass M is assumed to include the mass of the track and small mass m. The small 
masses are themselves connected by weaker springs. Assigning longitudinal displacements u„ to the large 
masses Af, and transverse displacements w„ to the small masses m in the frictionless tracks, the Hamiltonian 
of the mechanical model is: 



H = J2 l^ill + ^mvl + W{L + Un+l -Un) + V U{L + ul^^ - U„)2 + {Vn+1 - Vn)A 



(92) 



where W is the potential energy of the stiff springs connecting the large masses and V is the potential energy 
of the weaker springs. 

The associated equations of motion are 

Miln = W {L + Un+1 - Un) - W {L + Un - Un-l) 

+ S{Dn+l){L + Un+1 - Un) - S{Dn){L + Un - Un-l) (93) 
mVn = S{Dn+l){Vn+l - Vn) - S{Dn){Vn - Vn-l) , 

where we have set S{D) := V'{D)/D and Dn ■= \/{L + Un — Un-i)^ + (vn — Vn-i)^- It is clear that one 
may take the undisturbed state of the internal modes Vn+i = Vn for all n to hold exactly, in which ease only 
the backbone motion is relevant. We will be interested in small amplitude, linear motions of the Vn and how 
they are affected by the motion of the backbone. 

The disparity between the masses M and m, and that of the strengths of the associated springs, is 
introduced by letting /x be a small parameter and assuming m = fiM and V = nU (and correspondingly, 
S = fiZ). We make the small-amplitude long- wave ansatz: 

Unit) = hu{X = hn,T = ht) , and w„(t) = hv{X = hn,T = ht) , (94) 

where /i is a small lattice-spacing parameter. Expanding the functions W and Z in Taylor series about the 
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equilibrium position, the equations of motion become 



Mh^^d^u = h^W'\L)d\u+ —W'\L)dj^u + h^W''\L)dxu-d\u + 0{h'^)+0{h^lJL), 
Mh^d^v = h^Z{L)d\v + — Z(L)9^v + h^Z'{L)d\u ■ dxv + h^Z'{L)dxu ■ d\v + ©(/i^) . 



(95) 



Trapping of ?;-waves by w-waves becomes possible if the wave speeds are equal. Therefore, we assume that 
W"[L) = Z{L) = Mc^. Changing variables to x = ^ - cT and r = h'^T yields 



Mc 



Mh?dlu - 2Mcd^drU = -r^d^u + W"'{L)d,^u ■ d^u + 0{h) + 0{fi/h^) , 

12 



(96) 



Mh'^d^v-2Mcd^drV = -^d^v + Z\L)dlu ■ d^v + Z\L)d^u ■ d^v + 0{h) . 
As h [Q with n <^ h"^, we find the coupled system 



drA 



1 



—dlA 

24 ^ 



5v 



= 0, 
= 0, 



(97) 



as a formal limit, where A = W" {L)d^u/ {2Mc) and B = d^v. After a simple rcscaling of x t, this 
becomes (1) and (2) with k = Z'{L)/W"'{L). As described in § 6, the influence of solitons on linear waves 
can be qualitatively different for different values of the coupling constant k, with important bifurcations 
occm-ring at the values k = k}^^ = (n+ l)(n + 2)/12, for n = 1,2,3,.... As we have seen, this coupled system 
can be solved exactly in (at least) two cases: k = 1, and k = 1/2. The former case is just the linearized 
KdV; see Sachs [S83]. The latter case is the one that is solved in the main text of this paper. 
Consider this example with the potential of the strong and weak springs given respectively by 



W{D) := ^n^D' + ^aD^ , 



and ViD) := f^i^K^D^ + ^l3D^) ■ 



Thus (3 = 3aK and the condition that the wave speeds are equal is 



1 



aL^{K — 1) = Kw — Hv ' 



(98) 



(99) 



The effect of each bifurcation point in k is now clear. At k = 1/2 the first harmonic of the w- waves begins 
to resonate with the w-waves. As k increases through k = 1 we pass through a transition from supercritical 
resonance to subcritical resonance. Similarly, at the odd bifurcation points, k = K2to_i, for m = 1, 2, 3, . . ., 
the m*^ harmonic w-wave begins to resonate with the u- waves. Then at the even bifurcation points, k = Kj™, 
for TO = 1, 2, 3, . . ., the nature of the resonance for this mode changes from supercritical to subcritical. 

Coupled systems of equations like the pair (1) and (2) often arise as formal asymptotic reductions of 
mechanical models for complicated one-dimensional waves. Often these asymptotic models are integrable. 
For example, in an elastic rod, the interaction between axial twist waves and helical deformation waves gives 
rise to an integrable Manakov system of coupled nonlinear Schrodinger equations [LG99]. 

The coupled system (1) and (2) for k = 1/2 is also intimately connected with an integrable multicom- 
ponent (an arbitrary number of components, all appearing symmetrically) coupled KdV equation [MC99]. 
Indeed, from one point of view it is this connection that yields the solvability of (1) and (2) for k = 1/2 
described in detail in this paper. The solution method presented here also can be used to give the complete 
solution of the coupled KdV system. That system in turn can be interpreted as a phenomenological model 
for the transport of the mass integral through an A''-soliton solution of KdV [MC99]. 

Finally, we would like to point out that there are also some applications in which linear equations 
of the form (2) occur with KA{x,t) being a given function. In this case, c{x,t) = KA{x,t) represents 
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a given spatiotcnnporal modulation of the speed of linear dispersive waves, say due to propagation in an 
inhomogeneous medium. Such problems arise in the modeling of the propagation of weak internal waves in 
a channel of varying width [CG99]. For such applications, we may view the solvability of the coupled system 
(1) and (2) for k = 1 and k = 1/2 as a kind of (big) catalog of special cases of the function c(x,t) for which 
the linear equation (2) is solvable in its own right. For other values of k > 1/2 the linear wave system is 
unstable, while for all values of k < 1/2 it is stable. 
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